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Abstract 

The combination of collinear factorization with effective field theory originally developed for soft interactions of heavy 
quarks provides the foundations of the theory of exclusive and semi-inclusive B decays. In this article I summarize 
some of the later conceptual developments of the so-called QCD factorization approach that make use of soft-collinear 
effective theory. Then I discuss the status and results of the calculation of the hard-scattering functions at the next 
order, and review very briefly some of the phenomenology, covering aspects of charmless, electroweak penguin and 
radiative (semi-leptonic) decays. 
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1. Introduction 

In the early 1990s the heavy-quark effective theory 
ID El|3] and the heavy-quark expansion 000 were 
developed into powerful tools to describe the dynamics 
and emerging symmetries 0 [8| of QCD in the limit, 
when the mass nib of a quark is much larger than the 
intrinsic scale A of the strong interaction. The basic 
assumption for the validity of this expansion is that the 
limit is taken when all other dimensionful parameters, 
including the momenta of particles, are fixed to be of 
order A or another soft scale. After integrating out the 
scale nib, the effective dynamics is governed exclusively 
by the soft scale, since no other scale is present in the 
problem. 

In the second half of the 1990s the CFEO experiment 
at Cornell accumulated enough statistics of B mesons 
to measure for the first time the rare b —> u and loop- 
induced decays in exclusive two-body final states. The 
tools above do not apply to this situation. When the fi¬ 
nal state consists of two light mesons, say two pions, 
the energy of the pions is roughly nib 12. The presence 
of a large external momentum invalidates the assump¬ 
tion that all small-virtuality fluctuations are soft. Ener¬ 
getic, light particles can radiate energetic particles and 
remain close to their mass-shell when the radiation is 


collinear to their direction of flight. When the invariant 
mass squared is of order A 2 as must be the case for an 
exclusive final state of light mesons, collinear radiation 
cannot be described perturbatively. 

The idea of collinear factorization is central to the 
QCD treatment of high-energy scattering with classic 
applications to deep-inelastic scattering, exclusive form 
factors at high-momentum transfer, and jet physics. In 
these cases, a crucial step is often the demonstration 
that soft effects cancel, such that the non-perturbative, 
collinear physics can be isolated in parton distribu¬ 
tions and light-cone distribution amplitudes. The then 
novel experimental accessibility of exclusive, energetic 
final states in fi-meson decays required the extension 
of collinear factorization to situations where the large 
energy is injected into the process not by a colourless 
external source such as the electromagnetic current, but 
the weak decay of a heavy quark surrounded by the soft 
degrees of freedom that make up the B meson. Soft 
physics is therefore expected to be more relevant than 
in traditional applications of collinear factorization. 

The problem was first addressed in a systematic way 
by Buchalla, Neubert, Sachrajda and myself isnionni 
motivated by the search for a rigorous QCD treatment of 
exclusive heavy-light ( Dn ) and light-light ( nn, “charm- 
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Figure 1: Graphical representation of the factorization formula. Only 
one of the two form-factor terms in |TJ is shown for simplicity. Figure 
from llOl . 

less”) final states, based on nothing but the scale hier¬ 
archy mt, » A and the concepts of perturbative soft 
and collinear factorization. Assuming that the elec- 
troweak scale Mw is already integrated out, the start¬ 
ing point is the weak effective Lagrangian, which con¬ 
tains the flavour-changing interactions in the form of lo¬ 
cal dimension-six operators Q, of the four-fermion or 
chromo- and electromagnetic dipole type. The relevant 
objects to describe the B —> M\ Ms decay amplitude are 
then the matrix elements ( M\Mf\Qi\B ) of these opera¬ 
tors. For the further discussion below, I quote the fac¬ 
torization formula, valid in the heavy-quark limit up to 
corrections of order A/nib, in the form originally given 

m muni: 

{M x M 2 \Q\B) = 

X f 0 du T h (u) 

+ M 2 ) 

+ I dgdudv Tj l (aj, u , v) <p B (g) (f> Ml (v) 0 m 2 O) 

Jo 

if M\ and Ms are both light, (1) 

{M,M 2 \Qi\B) = 

£ du tJj(u) 4>m 2 (m) 

if Mi is heavy and Ms is light. (2) 

Here t ) denotes a B — * M\s form factor, and 

(/>x(u) is the light-cone distribution amplitude (LCDA) 
for the quark-antiquark Fock state of meson X. Tj.(u ) 
and 7j n (£, u, v) are hard-scattering functions, which are 
calculable perturbatively in the strong coupling a s \ m \2 
denote the light meson masses. Eq. 0 is represented 
graphically in Figure [T] The presence of a form fac¬ 
tor and of LCDAs shows that soft and collinear physics 
is relevant even at the leading order in the heavy-quark 
expansion. 

Eq- Q applies to decays into two light mesons, for 


which the spectator quark in the B meson can go to ei¬ 
ther of the final-state mesons. An example is the decay 
B —> n°K . If the spectator quark can go only to one 
of the final-state mesons, as for example in B c i —» n ' K , 
we call this meson M\ and the second form-factor term 
related to M\ «-» Ms on the right-hand side of 0 is 
absent. The factorization formula simplifies when the 
spectator quark goes to a heavy meson [see (|2j], such 
as in B,/ — > !)' n . In this case the hard interactions 
with the spectator quark represented by the third term in 
0 can be dropped, because they are power-suppressed 
in the heavy-quark limit. In the opposite situation that 
the spectator quark goes to a light meson but the other 
meson is heavy, factorization does not hold. 

The simplest case to gain some intuition is B,i —» 
D + n~, when the D meson is also assumed to be para¬ 
metrically heavy. The spectator quark and other light 
degrees of freedom in the B meson have to rearrange 
themselves only slightly to form a D meson together 
with the charm quark created in the weak b —> cud tran¬ 
sition. For the other two light quarks ud to form a pion 
with energy O(nib), they must be highly energetic and 
collinear, and in a colour-singlet configuration. Soft in¬ 
teractions decouple from such a configuration and this 
allows it to leave the decay region without interfering 
with the D meson formation tm The probability of 
such a special configuration to form a pion is described 
by the leading-twist pion LCDA <p n {u). The factoriza¬ 
tion formula |2| provides the quantitative framework for 
this discussion, which allows to compute higher-order 
corrections. For example, if the light quark-antiquark 
pair is initially formed in a colour-octet state, one can 
still show that soft gluons decouple at leading order in 
A /mb, if this pair is to end up as a pion. This implies 
that the pair must interact with a hard gluon, and hence 
this provides a calculable strong-interaction correction 
to the basic mechanism discussed above[j] An important 
element in demonstrating the suppression of soft inter¬ 
actions, except for those parametrized by the B —> D 
form factor, is the assumption that the pion LCDA van¬ 
ishes at least linearly as the longitudinal momentum 
fraction approaches the endpoints u = 0,1. This as¬ 
sumption can be justified by the fact that it is satisfied by 
the asymptotic distribution amplitude (f> n {u) = 6 u{ 1 - u), 
which is the appropriate one in the infinite heavy-quark 
mass limit. If this were not the case, the pion could be 
formed with larger probability in an asymmetric config¬ 
uration, in which one constituent has soft momentum of 
order A. Soft gluons would not decouple from this soft 


[ A correction of this type was computed already in fl3l , but the 
generality and importance of the result went unnoticed. 
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constituent, thus spoiling factorization. 

The above discussion relies crucially on the spectator 
quark in the B meson going to the heavy meson in the 
final state. If, as in the case of a D°7i° final state, the 
spectator quark must be picked up by the light meson, 
the amplitude is suppressed by the B —> n form factor. 
But since the D meson’s size is of order 1/A, the D° 
formation and B —» n transition cannot be assumed to 
not interfere and factorization is violated. 

The case of two light final state mesons is the most 
interesting one. The dominant decay process is indeed 
the same as for the case of the I ) 1 n final state, but this 
implies that the light meson that picks up the spectator 
quark is formed in a very asymmetric configuration in 
which the spectator quark carries a tiny fraction A/;;;/, 
of the total momentum of the light meson. Such a con¬ 
figuration is suppressed by the endpoint behaviour of 
the LCDA as discussed above, but owing to this sup¬ 
pression there exists a competing process, in which a 
hard gluon is exchanged with the spectator quark, pro¬ 
pelling it to large energy, thus avoiding the endpoint- 
suppression penalty factor. If the hard gluon connects 
to the quark-antiquark pair emanating from the weak 
decay vertex to form the other light meson, this gives 
rise to the third contribution to the factorization formula 
Q- This further contribution, called “hard-spectator in¬ 
teraction”, involves three LCDAs, including the one of 
the B meson and resembles the expressions that appear 
in the theory of exclusive form factors at large momen¬ 
tum transfer nuns. 

The above discussion reflects the theoretical under¬ 
standing of QCD factorization in exclusive B decays 
as of around 2000. The present article summarizes re¬ 
sults from a long-term project on soft-collinear factor¬ 
ization in B-me son decays started in 2003. The three 
following sections are devoted to three different strains 
of research. The first covers the further conceptual de¬ 
velopment of QCD factorization which centres around 
(a) factorization of the B —> n form factor and hard- 
spectator scattering, (b) a calculable example of factor¬ 
ization of endpoint divergences, and (c) factorization at 
sub-leading power in A/m/, for the semi-inclusive semi- 
leptonic decays B —> XJ.v. The hard-scattering kernels 
T\.(u) and 7j n ( u, v) were computed in @[I0) at 0{a s ), 
which corresponds to the one-loop and tree approxima¬ 
tion, respectively. Since then most of the required cal¬ 
culations for the next order have been done and will be 
discussed in the second part. One of the most remark¬ 
able implications of <[T}, ([2]) is that the strong phases, 
which must be present to make direct CP violation ob¬ 
servable, are either 0(a s ) from loops at leading power, 
or 0(A/nib) power-suppressed. The next order there¬ 


fore also yields the first QCD correction to the direct CP 
asymmetries, which is usually required to obtain a reli¬ 
able result. Finally, some aspects of the phenomenology 
of the factorization approach for charmless and other 
exclusive decays will be very briefly presented in the 
third section. 

I conclude this introduction with an apology. The cir¬ 
cumstances that require the writing of this article imply 
that it focuses on my own work. It does not do justice 
to that of many others and does not substitute a review 
of the subject that remains yet to be written. 

2. Conceptual development of QCD factorization 

The further conceptual development of QCD fac¬ 
torization in B decays has greatly benefited from the 
development of soft-collinear effective theory (SCET) 
maim mi Era, the effective Lagrangian formula¬ 
tion of diagrammatic factorization. In turn, the desire to 
better understand the factorization of exclusive heavy- 
quark decays has been instrumental in the early develop¬ 
ment of SCET. For example, the essence of the diagram¬ 
matic two-loop analysis IfTOl that demonstrated that the 
structure of infrared singularities in the quark represen¬ 
tation of the B —> Dn amplitude is consistent with <§, 
is elegantly reproduced by the decoupling of soft glu¬ 
ons from the leading-power collinear SCET Lagrangian 
through a field redefinition by a soft Wilson line, and 
the subsequent cancellation of the Wilson lines from 
the colour-singlet current that overlaps with the pion 
state ED- This proves factorization to all orders, pro¬ 
vided that SCET reproduces all the relevant infrared de¬ 
grees of freedom, which is the standard assumption. In 
the following, however, I focus on hard-spectator scat¬ 
tering, where not all the factorization properties are evi¬ 
dent in (JT|), since two perturbative scales nib and fin/, A 
are hidden in the scattering kernel T?(£, u, v). It turns 
out that the main issues can already be understood from 
the factorization of heavy-to-light form factors, so I turn 
to these first. 

2.1. Heavy-to-light form factors 

The heavy-to-light form factors parametrize matrix 
elements of the form {M(p')\qTb\B(p)), where T de¬ 
notes a Dirac matrix, and we are interested here in the 
region of large momentum transfer of 0(mb ) to the light 
meson M. The physics contained in this matrix element 
is surprisingly rich and complex. First attempts com¬ 
pute the large-recoil form factor in terms of LCDAs date 
back to l22l and are based on the assumption that the 
light meson is produced in a configuration where the 



4 


M. Beneke, Soft-collinear factorization in B decays (2015) / |2 7 | 


s 


s 


he 


c 


c 



Figure 2: Tree diagrams for the B —> n form factor. Left and middle: 
the pion is produced in a symmetric configuration. In the middle dia¬ 
gram hard fluctuations have already been integrated out. Right: when 
the pion is produced in an asymmetric configuration, the diagram is 
purely soft. 


Table 1: Scaling of momentum modes relevant to B —> M form fac¬ 
tors. 


{n+p,pi_,n-p) 

terminology 

virtuality 

(U,l) 

hard 

1 

(1, d, A 2 ) 

hard-collinear 

A 2 

(A 2 , A 2 , A 2 ) 

soft 

A 4 

(1 ,d 2 ,d 4 ) 

collinear 

A 4 


quark and antiquark constituents both carry large mo¬ 
mentum, which requires the exchange of an energetic 
gluon with large virtuality m b A (left diagram in Fig¬ 
ure [2}. However, the calculation exhibits an endpoint 
divergence, which invalidates this assumption. The op¬ 
posite starting point was adopted in (23l . where it was 
shown that the three (seven) independent scalar form 
factors that parametrize the transition of a pseudo-scalar 
B meson to a pseudo-scalar (vector) light meson can be 
expressed in terms of only a single (of only two) func- 
tion(s), if the transition is dominated by soft interac¬ 
tions. This corresponds to the right diagram in Figure [2] 
and implies that the light meson is always produced in 
an asymmetric configuration, which though improbable 
is favoured by the presence of the soft remnant of the 
/1-meson after the heavy-quark decay. Ref. Il24ll showed 
that the symmetry relations that leads to the reduction of 
the number of independent form factors are broken by 
perturbatively calculable corrections, and that the end¬ 
point divergences that appeared in the hard-scattering 
approach could be absorbed into the soft form factors. 
This led to the conjecture that at leading power in the 
A /nib expansion the scalar form factors factorize as 

F* M (q 2 ) = Cfq 2 ) £ BM (q 2 ) +cf B ® T,(q 2 ) ® <f> M (3) 

at large recoil (24l . Here <S> denotes a convolution in 
momentum fraction, and %bm(t) refers to the soft form 
factor that satisfies the symmetry relations. 

With the advent of SCET the heavy-to-light form fac¬ 


tors were extensively studied H3[m|22|2S]|27||2!| re¬ 
sulting in a “proof” of |3]l in l27l . We first note that 
there are two large scales in the problem. Their ex¬ 
istence can already be seen in the left tree diagram in 
Figure [2] In order to convert the soft spectator quark 
with momentum A into a energetic light quark, the ex¬ 
changed gluon must have virtuality itibA. The heavy- 
quark line that connects the gluon vertex to the external 
weak vertex, however, is off-shell by the larger amount 
m 2 b . Factorization should therefore proceed in two steps 
through the construction of two effective theories, above 
and below the scale yfmbK (26). When the off-shell 
heavy quark in the left diagram is integrated out, the 
effective interaction looks as in the middle diagram of 
Figure [2] which suggests that higher-dimensional oper¬ 
ators with additional gluon fields must be relevant (26). 

Let the four-momentum of the light meson be p' = 
Eti-, where n± define two light-like momenta with 
n+ ■ n- = 2, and A = yjA/mt, defines the small pa¬ 
rameter that organizes the power-counting in SCET The 
relevant momentum modes are summarized in Table Q] 
In the B-meson rest frame, the B meson consists of a 
heavy quark with soft residual momentum and soft con¬ 
stituents, while the light meson is made of collinear par- 
tons. The virtuality of soft and collinear modes is A 2 
and non-perturbative. In contrast, the hard and hard- 
collinear modes can be integrated out perturbatively one 
after the other. Each step defines a new effective theory, 
called SCETi and SCETn in (26), and SCET(hc,c,s) and 
SCET(c.s) in 127) . In the following I summarize the key 
results in each step, referring to the original papers for 
many important details. To be definite, M will be as¬ 
sumed to be a pion. 

In the first step, the hard modes are integrated out. 
The SCET Lagrangian is not renormalized to any order 
ed, but the QCD currents <//I Q must be matched to 
SCETi. The relevant SCETi current operators have field 
content (26) 

%cK, (4) 

where here c can be c or he. The two operators have 
different A scaling in SCETi, but nevertheless both con¬ 
tribute to the leading power in A /mb to the form factor. 
The reason for this is that SCETi contains modes of dif¬ 
ferent virtualities, so that the scaling of the matrix ele¬ 
ments of ([4ji is determined only after an analysis of the 
time-ordered products of the currents with interactions 
from the SCETi Lagrangian. Power-suppressed inter¬ 
actions are needed to obtain a term with non-vanishing 
overlap with the pion and /1-meson state, and the analy¬ 
sis shows that both types of operators contribute at order 
( A/m b ) 2 ' 2 ( 23 . 
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Since the (hard-) collinear modes have n+-p ~ 0(m b ), 
SCET operators are generally non-local and have to be 


dressed with collinear Wilson lines W c to make them 
gauge-invariant. The precise matching equation is 


(0T,■ 0(0) = J ds C\ A0 \s) O (A0 \s ; 0) 


J dhdh 1, S 2 ; 0) + • • •, 


(5) 


where 

O (A0) ( S ; x) = (£W c )(x + sn + )h v (x _) = QW c ) s h v , 0 B1) ( Sl , s 2 ; x) = — (fW c ),, (Wj iD ±c/J W c ) S2 h v (6) 

m b 

correspond to the two operators in (Jdj), and xf‘_ - (n+ ■x)n l ‘ /2, si = Sj>n b /n v. One of the ^-integrations can be removed 
upon taking the matrix elements by using translation invariance. We then define two leading-power SCETi pion form 
factors through 

(n(p'MW c )h v \B v ) = 2 £?*,(£), (7) 


W)\ —QWcXW] iip cl _ W c )(rn + )h v \B v ) -2 E f 
m b J 


dre 


ilErr -7 


E Bn {T,E). 


( 8 ) 


Here |B,,) denotes the B-meson state in the static limit normalized to 2»ij. At this step, the three pion form factors can 
be represented as 


Ff*(E) = Ci(E) £ B7r (E) 


f 


dTC ( i m> (E,r)E B AT,E), 


i = +, 0, T. 


(9) 


The point to note here is that the three form factors of 
the B —» n transition can be expressed in terms of a sin¬ 
gle form factor £n n (E) and another non-local form factor 
S b „(t,E). A number of “symmetry” relations between 
form factors emerges already at this stage. 

The first term of (|9]( is already in the form of Q and 
is often associated with the “soft form factor”. While it 
certainly includes the soft overlap contribution, i.e. the 
contribution where the pion is formed in the asymmet¬ 
ric configuration (see the right diagram of Figure [2}, it 
is defined as a SCET[ matrix elements still containing 
hard-collinear short-distance effects. The matching of 
the operator QW c ) s h v to SCETn by integrating out the 
hard-collinear effects has been studied in [27. 281. but 
has so far been unsuccessful. The problem arises from 
the absence of a consistent interpretation of endpoint di¬ 
vergences that appear in the convolution integrals of the 


matching coefficients, which point to missing degrees of 
freedom. The analysis of j27l shows that the structure 
must be complicated, since three-particle LCDAs of the 
/i-meson and pion appear at leading power. £ Bj! (E) is 
therefore treated as an unknown, non-perturbative func¬ 
tion. Fortunately, there is only one instead of the origi¬ 
nal three (two instead of seven for vector mesons). 

Yet (|9| would be of little use, if the second term con¬ 
taining an unknown form factor E Bji (t, E) of two vari¬ 
ables could not be simplified. In the second match¬ 
ing step, one therefore integrates out the hard-collinear 
scale y/m b A. Note that this step is only applied to the 
second term and not the first for the reasons mentioned 
above. A power-counting analysis EH shows that the 
operator O 1 ^ 1 ’ ~ % c A^ hc h v matches only on four-quark 
operators of the form (q s h v )(^ c ^ c ) with no additional 
fields or derivatives at leading power in the A expansion. 
The explicit operator matching relation reads 


-/ 


2 E I — e~ l2ETr 

2n 


(i WcXOXWj i#) c± W c )(rn + )h v ( 0) 


dcodv J(t; v, ln( Eu>/p~)) 


(tW c )(sn + )^y 5 (W^m 


Jft 


(.q s Y s )(tn-)’i-y 5 (Ylh v )( 0) 

2 ‘ Jft 


( 10 ) 










with 
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G(V) 

/ J (w) 


(fw c )(5/i + )^r 5 (wj^)(0) 

{qJ s )(tn-)l—y 5 (Y' s h v )(()) 


^ J ds e~ ,svn + p ' (fW c )(sn + )^-y 5 (W^m, 
Y n f dte i,u (q s Y s )(tnJ-y 5 (Y%m, 


(ID 

( 12 ) 


and 7(r; v, \n(E(L>//i 2 )) the hard-collinear matching co¬ 
efficient, which can be computed perturbatively. (The 
middle diagram in Figure [2] is one of the two tree di¬ 
agrams that contribute to J.) Y s denotes a soft Wil¬ 
son line. The Wilson lines W c , Y s ensure the gauge- 
invariance of the light-cone operators. 

In the SCET|i Lagrangian there are no interactions 
between collinear and soft fields, since the sum of a 
collinear and soft momentum has hard-collinear virtu¬ 
ality, which is already integrated out. The SCETn ma¬ 
trix element of a product of collinear and soft fields as it 
appears on the right-hand side of ([TO]) always falls apart 


into separate collinear and soft matrix elements. In the 
present case, these define the LCDA 

W)lg(v)IO) = -if n E Mv) (13) 

for the pion and 

<0|PM|5 V > = <p B+ (oj) (14) 

for the /f-meson MM- Thus the four-fermion op¬ 
erator factorizes into simpler matrix elements upon in¬ 
tegrating out the hard-collinear scale and the non-local 
form factor can be written as m 


dv J(t\ v, MEoj/fi 2 )) f B <l>B+(oS)M„(y). (15) 


- / I7\ lllB 

*Bn(T, E ) = -- 

4 mt, 


s-»oo /-»1 

( dto 

Jo Jo 


An important point needs to be made here. Due to the 
absence of soft-collinear interactions in the SCET n La¬ 
grangian any operator at any order in the A expansion 
formally factorizes into a convolution of a soft and a 
collinear function, but the result is not always correct, 
since the convolution integrals often diverge at the end¬ 
points — see the discussion of ^Bn(E) above. The fi¬ 
nal part of the factorization “proof” therefore consists of 
showing that the convolution integrals in ( p~5j ) converge. 
The argument given in l27l shows that J can depend on 
(o only through terms of the form 1/w x In "(Eaj/fj 2 ), 
which together with the endpoint behaviour of the B- 
meson LCDA guarantees the convergence of the co- 
integral. From this the convergence of the v-integral is 
inferred, since an endpoint divergence in a collinear in¬ 
tegral must always have a correspondence in a soft con¬ 
volution. 

The subtleties of soft-collinear factorization are re¬ 
lated to the factorization of modes with equal rather 
than hierarchical virtuality ( 1271 . Figure 2), which are 
distinguished by the scaling of a longitudinal light-cone 
momentum component rather than transverse momen¬ 
tum. The factorization of such momentum regions usu¬ 


ally leads to divergences that are not regulated in dimen¬ 
sional regularization. Other schemes that work such as 
analytic regularization whereby one raises certain prop¬ 
agators to non-integer powers l/[-£ 2 ] 1+fl violate some 
of the longitudinal boost symmetries of the classical La¬ 
grangian, which constrain the dependence of coefficient 
functions and operators on light-cone momentum com¬ 
ponents in the naive factorization formula. The absence 
of a regulator that preserves all the properties of the clas¬ 
sical theory at the quantum level implies that SCETn 
factorization is generally “anomalous’ ’ HO) and valid 
only in special cases, such as the form factor E Brr (r,E) 
discussed above. To my knowledge, a general argument 
when SCETn factorization is not anomalous is not avail¬ 
able. These difficulties do not affect the first matching 
step where hard modes are integrated out. At least, I do 
not know of an example for which the convolution in 
the hard-collinear momentum fraction (r above) is di¬ 
vergent or dimensional regularization fails in matching 
the full theory to SCETn 

Returning to the heavy-to-light form factor we plug 
© into (|9|) and obtain the conjectured factorization 
formula |3|, provided we identify the kernel T,(q 2 ) of 
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r ! (£;v,lnw)= -p- I dTCf l, (E,T)f B f n J(T-vM(Eajhi 1 )), 
4rrib 


(16) 


or Ti = C tm 1 ® 7 in short. The hard-scattering ker¬ 
nels of spectator scattering are themselves a convolution 
of a hard and a hard-collinear function l26l . The one- 
loop corrections to C\ Bl) 113 1 1 [32] and J lf33l l34l have 
both been calculated, and the explicit results confirm the 
general factorization formula, including the arguments 
supporting the convergence of the convolution integrals. 
We therefore conclude that the three independent B —> n 
form factors can be expressed in terms of a single form 
factor g Bjr plus corrections that begin at 0(a s ( sJm h A)) 
and can be expressed in terms of the universal, process- 
independent LCDAs of the mesons. 

2.2. Spectator scattering in hadronic decays 

I have discussed the heavy-to-light form factors at 
length, since after these preparations the extension to 
charmless B decays is relatively simple. Charmless de¬ 
cays were first considered in SCET in H351 [36). After 
the development of the theory of spectator scattering 


for the form-factor, the derivation of the QCD factor¬ 
ization formula within the SCET framework was com¬ 
pleted in f37j . 

Instead of current operators one needs to match the 
operators present in the weak effective Hamiltonian. As 
an example, consider the contribution of Q - [TihYi 1 - 
y^baMay^ -75 )ub\ to the decay B d -> tx + n~ (a.b are 
colour indices). The extra n relative to the B —* n tran¬ 
sition has large momentum q in the direction opposite 
to 7r + . The SCET Lagrangian is therefore extended by a 
second collinear sector with fields x- ^c 2 in addition to 

A c i, satisfying = 0 and t/+x = 0. Analogous to 

0 the operator Q matches to the expression 


Q= JdtT 1 (t)0 1 (t)+ JdtdsH ll (t,s)O ll (t,s) (17) 


at leading power after integrating out the hard modes, 
where the two operator structures are now given by 


0\t) = (xW C 2 )(tn-y—(\ - y 5 )(Wl 2 x) 


i-7s)Av 


f dsC'f I, (s)(^W cl )fi + [W^ 1 ip ±cl W cl ](sn + )( 1 +y 5 )h v 
m b J 7+ 


MB 1), 


O l \t,s) = — 


1 

m b 


(xW c2 )(tn-) — ( 1 -y 5 )(W' c2 x) 


(IWcO^twJ^xdWdKsn+X 1 + y s )hy 


( 18 ) 

(19) 


It might appear that a final state with two pions should 
create a major problem, since the physics of charmless 
decays is obviously more complicated than that of the 
form factor. However, the two collinear sectors are al¬ 
ready decoupled after the first matching step to SCETi, 
since the sum of a cl and c2 momentum necessarily 
has hard virtuality 0(mfy. Thus there can be no cou¬ 
pling of the two modes in the SCET] Lagrangian, and 
hence the (7r~7r + |... \B) matrix element of the opera¬ 
tors ( p~8] >, ( [T9| factorizes already at the hard-collinear 
scale into a (7r~|... |0) matrix element of the c2 fields 
and the (7r + |... \B) matrix element of the remainder af¬ 
ter decoupling the soft gluons from the c2 fields by the 


soft Wilson-line field redefinition in the same way as for 
B d —> D' ti . The former simply gives the pion LCDA 
(f>M 2 (u) in 0 . The latter matrix elements are the same as 
those that appear in the factorization of the B —> n form 
factors. Following f38l . the square bracket in ( | 1 8[ > has 
been defined such that its matrix element coincides with 
the full QCD B —> n form factor. This yields the first 
term in 0- The matrix element of the second square 
bracket in ( [T9| can be related to 2L Bn defined in <|8), 
so ( p~5] > can also be used here. Introducing the Fourier 
transform 

H n (u, v) = f dids e ,<ui+( '- m 


( 20 ) 
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of the hard-matching coefficient of the operator O u (t, s ) respond to collinear momentum fractions, one finds that 

associated with the spectator scattering, where u, v cor- the kernel T u (a>, u, v) in (jTJ is given by the convolution 


r* 1 

T 11 (oj, u, v) = — I dzH u (u,z) J(l -z;v, of). 

8 m b Jo 

Here J is the same hard-collinear matching coefficient that appeared in the factorization of the form factor, 
reproduces the QCD factorization formula (JT]» for charmless decays in the SCET framework. 


( 21 ) 

This 


The SCET derivation yields several new insights. 
First, the spectator-scattering kernel is a convolution 
7’ 11 = H 11 <8> J of a hard and hard-collinear function, 
and precise operator definitions can be given, which are 
important for consistent higher-order calculations (see 
following section). Second, the decoupling of the sec¬ 
ond light meson at the hard-collinear scale implies that 
the perturbative strong interaction phases originate from 
the imaginary part of hard loops, and never from hard- 
collinear loops. Third, the upper limit of the ^-integral 
in {T} should be replaced by an integral over to with up¬ 
per limit oo instead of 1. The variable J = to/m B corre¬ 
sponds to a light-cone momentum fraction of the spec¬ 
tator quark in the B meson. However, the LCDA <p B (co) 
is defined in SCET n , which contains a static quark field 
h v of formally infinite mass. Relative to this the variable 
io can take any value. The one-loop renormalization of 
tf> B (to) ll39l shows that a perturbative tail is generated 
that ranges to to = oo, unless the LCDA is defined in a 
cut-off scheme. Since the characteristic scale is to ~ A, 
the B-meson LDCA can appear in leading-power fac¬ 
torization formulae only in the form of the inverse mo¬ 
ments 

J r>o° j 

— (p B +(co,ld) In” —, (22) 

o oj to 

where the most important moment A B (/o) is defined by 
<xo(/u) = 1 . At order 0(a k s ), up to n - 2k - 2 logarithms 
can appear, so the number of parameters related to the 
B-meson LCDA is limited to a few in practical calcula¬ 
tions. 

2.3. Variants of factorization for hadronic decays 
Several variations of factorization as described above 
have been advocated. The PQCD framework ESI 
ED assumes that the B-meson transition form factors 
F bm '( 0) are also dominated by short-distance physics 
and factorize into light-cone distribution amplitudes. 
Both terms in (|T]» can then be combined to 

(MiM 2 \Qi\B) = T PQCD +4> b + <f> Ml * (23) 


The implicit assumption that the heavy-to-light form 
factors do not receive any soft long-distance contribu¬ 
tion (due to Sudakov suppression) seems difficult to jus¬ 
tify to me, and certainly contradicts the picture devel¬ 
oped above. At the technical level, the endpoint diver¬ 
gences are regularized by keeping the intrinsic trans¬ 
verse momentum of the meson constituents. There is 
a larger sensitivity to perturbative corrections at low 
scales, where the strong coupling is large and per¬ 
turbation theory potentially unreliable. From a phe¬ 
nomenological perspective, PQCD needs fewer non- 
perturbative input parameters, but there is a larger de¬ 
pendence on unknown light-cone distribution ampli¬ 
tudes. A principal difference between the PQCD and the 
QCD factorization approach described above is the rel¬ 
ative importance of the weak-annihilation mechanism 
and the generation of strong rescattering phases. In 
QCD factorization the phases arise at the scale m B from 
loop diagrams that have yet to be included in the PQCD 
approach, and from a model for power-suppressed weak 
annihilation, whereas in the most widely used imple¬ 
mentation of PQCD the strong phases originate only 
from a weak-annihilation tree diagram. A recent at¬ 
tempt to compute radiative corrections in the PQCD 
approach resulted in uncancelled infrared divergences 
m, which led to the introduction of final-state spe¬ 
cific non-perturbative parameters, in violation of factor¬ 
ization. Since a clear power-counting scheme in A/nib 
has never been established for the PQCD formula ( |23) >, 
this is not necessarily in contradiction with factoriza¬ 
tion at leading power in A/»;*. However, it prevents a 
systematic improvement of PQCD calculations beyond 
tree-level. 

Ref. 03 proposed a phenomenological implementa¬ 
tion of factorizatior0 that differs in two important re- 


2 Sometimes called the “SCET approach”, which I find mislead¬ 
ing, since the theoretical basis of diagrammatic QCD factorization 
and SCET factorization is exactly the same. See the derivation of (JTJ 
in the previous subsection. 
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spects from the BBNS approach EEBUU. First, 
perturbation theory at the intermediate scale 5 /nib A is 
avoided by not factorizing the spectator-scattering term 
into a hard and jet function. Eq. G then takes the form 

{M x M 2 \Qi\B) = f bm ' Tj * (f> M2 

+ E hm ' * H\' * d>M 2 , (24) 

which is simply the matrix element of the SCETj 
representation ( [17] ) of the four-quark operator of the 
weak effective Hamiltonian. Second, penguin diagrams 
with charm loops, the so-called “charming penguins” 
Il43l , are claimed to be non-factorizable, hence non- 
perturbative at leading order in the A/nib expansion. 
From the phenomenological perspective the principal 
difference to the BBNS approach concerns again the 
generation of strong interaction phases. Since the non¬ 
local form factor is unknown, Eq. ( |24| ) can be used in 
practice only at tree level, hence the matrix elements 
have no imaginary parts. The only exception is the ma¬ 
trix element corresponding to the charm penguin am¬ 
plitude, which is considered as an unknown complex 
number (one for each final state, to be fitted to data), 
and therefore represents the only source of direct CP vi¬ 
olation. Weak-annihilation effects are neglected, since 
they are power-suppressed. A critique of this approach 
is given in SU- From the theoretical perspective of fac¬ 
torization, the most important issue raised is the ques¬ 
tion whether the penguin loops with charm factorize or 
not. The concern is that for massive quark loops the 
non-relativistic energy scale m c v 2 near threshold is of 
order A, and hence there appears to be a sensitivity to 
non-perturbative scales not present for massless quark 
loops. The suppression of sensitivity to this region and 
the associated issue of quark-hadron duality for charm- 
quark loops has been discussed in P3ll . which resolves 
the issue in favour of factorization. 


2.4. Endpoint divergences in quarkonium final states 

Two-body decays to a charmonium (M 2 ) and a light 
meson (M\) were argued IfTUl to also factorize ac¬ 
cording to G in the non-relativistic (Coulombic) limit 
m c v 2 » A, since the size of the charmonium 1 /(m c v) 
is small compared to the wave-length of soft gluons, 
which therefore decouple. However, when the cal¬ 
culation of the hard-scattering functions was done for 
P-wave charmonia Bbl . uncancelled divergences ap¬ 
peared, including an endpoint divergence in spectator 
scattering, which contradicted the expectation of factor¬ 
ization. In this section I briefly review the result of P71 . 
which explains the origin of the divergences and resur¬ 
rects factorization at leading order in A/mt,. The main 
conceptual interest in this result derives from the fact 
that (to my knowledge) it constitutes the only example 
of a consistent factorization of an endpoint divergence 
in terms of precisely defined operator matrix elements. 

An important property of final states with charmo¬ 
nium when the latter is assumed to be a non-relativistic 
bound state is that there exist two sets of low-energy 
scales, m c v,m c v 2 related to the non-relativistic expan¬ 
sion, and y/nibA. A, related to the collinear expansion 
and the strong-interaction scale. In the following we as¬ 
sume that m c » y/mt,A » m c v, m c v 2 » A. 

Integrating out the scales nib, m c leads from QCD 
to an effective theory that combines SCE'lj with the 
non-relativistic effective theory of QCD for the charm 
quarks. A QCD operator such as Q — [cby'f 1 - 
y 5 )ba][s a 'y l j('\ - 75 ) 0 ,], which contributes to the decay 
Bit —t XcjK °, matches to SCE'lj® NRQCD operators 
which have exactly the same form as <[T8]), except 
that the part referring to the second light meson has to 
be replaced by a non-relativistic P-wave colour-singlet 
rather than a light-cone operator: 


(xW C 2)(tn-) y(l - 75)(V4C) 


0( 2S+1 P^) = <//% 



D t Av 


(25) 


with D 1 /- = I) 1 ' — (v • D)v M and E /( a spin matrix. One now 
finds that the tree-level matching coefficient of the oper¬ 
ator 0 11 (19 1 contains terms of the form m 2 l(m 2 b (\—y)) at 
tree-level, which are absent for massless quarks. Since 
y = 1 -y is the fraction of hard-collinear momentum car¬ 
ried by the gluon field in (19 1 , the divergence at y = 1 
occurs in the region when no momentum is transferred 


to the spectator quark, which is the endpoint region. In¬ 
deed, one then finds that the hard-spectator amplitude 
contains an integral of the form 


3\ 


hard spectator 
B—>H( ls+l Pj)K 


D 


.khms f 1 (pg(y) 

m b A B Jo ^ 


( 26 ) 


which diverges at y = 1 . 
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A crucial point is that the presence of the non- 
relativistic scales allows additional operators to con¬ 
tribute at leading order in A /nib. New and central to the 
present discussion are colour-octet operators such as 

0'j 3 Sf ) )=£W c y^(] -y 5 )T A h w ^T A Xv . (27) 

In the case of charmless decays to two light mesons the 
matrix elements of colour-octet operators can be non¬ 
zero only due to power-suppressed soft-gluon interac¬ 
tions, where soft means momentum of order A, thus 
they can be neglected at leading order in the A/nib ex¬ 
pansion. For charmonium, however, the decoupling of 
gluons with small momentum holds only when the mo¬ 
mentum is much smaller than m c v 2 ; gluons with mo¬ 
mentum m c v 2 contribute to the octet operator matrix el¬ 
ements even at leading order in A/nib. These contribu¬ 
tions are sub-leading in v, but so are the P-wave opera¬ 
tors due to the extra derivative in 0( 2S ' 1 Py ’). hence the 


gluon-exchange contribution to the S-wave octet opera¬ 
tors is relevant at leading order in the velocity expansion 
to P-wave charmonium production. 

Spectator scattering thus consists of hard spectator 
scattering, where the gluon exchanged between the cc 
system and the spectator quark in the B meson has en¬ 
ergy of order mb, and y ~ <9(1); and soft spectator 
scattering with gluon energy 0(m c v 2 ), in which case 
y ~ 0 (v 2 ) <s 1 is naturally in the endpoint region. 
Soft spectator scattering leaves the charm quarks close 
to their mass-shell and therefore does not reduce to 
an effective interaction of the form ( |~i~9| ). It is part of 
the colour-octet matrix element to be computed within 
the effective theory and is sensitive to the charmonium 
bound-state dynamics. The explicit calculation in the 
Coulombic limit m c v 2 » A gives an expression that 
contains 


soft spectator Ik/bHIb 
mAB 




dy<Pf;(y) 


7b 


2 y/zEh 
m b ( 1 -z)j m b Vl - Z 


(28) 


where E H < 0 is the binding energy of the charmo¬ 
nium state, z = 4m 2 /m 2 and y B = (m c a S C f) /2 is the 
inverse Bohr radius. We now compare the integral over 
the kaon distribution amplitude to ( |26} . While the in¬ 
tegrand there was applicable to y not near 1 and exhib¬ 
ited a logarithmic endpoint divergence as y —> 1, the 
integrand of (28 i is appropriate only to 1 — y ~ v 2 , i.e. 
in the endpoint region. There is no divergence here as 
y —> 1. However, for y » v 2 the integrand has the same 
logarithmic behaviour f dy(/>g(y)/y 2 as does the hard- 
spectator contribution for y <sc 1. Eq. (26 1 is based on 
approximations that are invalid when y is small, which 




soft spectator 
B—*H( 2s+l Pj)K 


D 


Zk/bHIb 

m b A B 


<pg( 1)( - ln/r - ln(l 


where A F is a real number given in Il47l depending on 
nib, ln c, Eh an d 7b- The endpoint contribution is pro¬ 
portional to <pf( 1), the derivative of the kaon LCDA at 
the endpoint, because the endpoint region is of size v 2 
rather than A /nib, hence it is justified to describe the 
quarks in the kaon by collinear quark fields. 


causes the endpoint divergence. We can regulate the di¬ 
vergent integral by cutting off the y integral above 1-fi. 
This corresponds to a hard factorization scale in the en¬ 
ergy of the gluon that connects to the spectator quark. 
The spectator-scattering contribution to the colour-octet 
matrix element originates precisely from the energy re¬ 
gion that is then cut out in ( |26| , and is not valid when 
y ~ (9(1), thus the correct interpretation of the y-integral 
in (28 i is J () dy —> J ay. To combine with (261 we 
must evaluate the regularized version of 
of order We then find 


up to terms 


z) + 1 + in + AT 7 ), 


(29) 


There are several remarkable features of this result. 


First, although (281 comes from a tree-diagram contri¬ 
bution to the octet matrix element, it carries a sizable 
absorptive part related to rescattering, which originates 
from the on-shell charm propagator in the integral over 
the charmonium bound-state wave function. Second, 
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the factorization scale dependence on /j cancels exactly 
with a corresponding logarithm in the (regulated) hard 
spectator contribution ( p6| ). The cancellation occurs as 
usual between an infrared (endpoint) divergence from 
the coefficient function and an ultraviolet (large energy) 
divergence of the octet matrix element. It further occurs 
between objects that have precise operator definitions in 
the relevant effective field theory. 

Unfortunately, these results do not readily generalize 
to final states of two light mesons. In the charmonium 
case above, there are two endpoint regions, y ~ v 2 and 
y ~ A/ttib, and we showed the cancellation of the end¬ 
point divergence in the former, while still working at 
leading order in A /mb- One difficulty with endpoint fac¬ 
torization at y ~ A/nib is that a light meson in the asym¬ 
metric qq configuration cannot be made out of generic 
collinear and soft modes in SCETn, since the sum of 
collinear and soft momenta has virtuality nib A, not A 2 . 
It seems that SCETn is lacking some of the relevant de¬ 
grees of freedom to write down operators that can de¬ 
scribe the endpoint region and overlap with the light- 
meson wave function. The understanding of this issue 
seems central to me to further theoretical progress on 
exclusive B decays, which has to address the factor¬ 
ization properties of power corrections in A /nib. The 
present use of factorization in exclusive B decays com¬ 
bines a rigorous theory at leading power in A /nib with 
a rudimentary parametrization of some power correc¬ 
tions, which introduces a considerable but presently un¬ 
avoidable arbitrariness in some quantities, especially di¬ 
rect CP asymmetries, where power-suppressed phases 
compete with loop-induced ones. 

2.5. Power corrections to semi-inclusive B decay 

While a rigorous theory of power-suppressed effects 
is not available for exclusive B decays, there is one im¬ 
portant class of processes for which soft-collinear fac¬ 
torization has been worked out completely including the 
0(A/mb) terms. These are the inclusive semi-leptonic 
decays B —> X u (v. where X u denotes a hadronic final 
state without charm, and inclusive radiative and elec- 
troweak penguin decays B — > X s y, B —> X s { + {~ in 
the kinematic region where the hadronic final-state is 
collimated into a single jet, which carries large energy 
of 0(nib), but has small (though not too small) invari¬ 
ant mass nib A ^ rn^. Leading-power factorization into 
hard, jet and soft (here called shape) functions for these 
processes was discussed in IfTTl |48| in the SCET and 
diagrammatic framework. 

The reason why a rigorous treatment of power cor¬ 
rections is possible here is the absence of collinear 
physics at the scale A. The two-step matching at the 


hard and hard-collinear scale then corresponds to the 
sequence QCD —> SCETi —> HQET of effective theo¬ 
ries. The jet functions containing the collinear physics 
at the scale yJnibA can be computed perturbatively in 
the strong coupling, while the non-perturbative soft 
physics is contained in well-defined HQET matrix el¬ 
ements. There is no other mode with virtuality A 2 
than the soft mode, when the virtuality of the inclu¬ 
sive hadronic final state is <9(ni/,A). The effective-theory 
framework provides a transparent book-keeping of the 
relevant interactions at every step in the calculation, in¬ 
cluding power-suppressed effects, and allows us to iden¬ 
tify all the relevant HQET operators. The factorization 
of semi-inclusive decays at 0(A/mb) has been discussed 
in Ii49ii50 . 511 at various levels of detail. The discussion 
below outlines the main points following the presenta¬ 
tion in ED, which contains further details and the ex¬ 
plicit tree-level computation. 

The strong interaction effects in semi-leptonic B de¬ 
cays are contained in the hadronic tensor VV 7 ' 1 ', related 
by the optical theorem to the imaginary part of the for¬ 
ward scattering amplitude, which we define as 

= - Im <B(v)|rHB(v)>. (30) 

n 

The correlator 7 7 ' 1 ' is the time-ordered product of two 
flavour-changing weak currents 

<3» 

where q is the momentum carried by the outgoing lepton 
pair in the B —» XJ’v decay and J 1 ' — qy 1 '(1 - y^)h. 

The matching of the current J^ from QCD to SCETi 
is the same as was discussed for exclusive decays. Due 
to rotational invariance in the transverse plane, a frame 
can be chosen such that the hard-collinear transverse 
momenta of order A = VA / m h appear only in the inter¬ 
nal integrations over the momenta in the jet, and since 
the integral over an odd number of transverse momenta 
either vanishes or must be proportional to one of the ex¬ 
ternal soft transverse momenta of order A 2 , the resultant 
expansion is in powers of A 2 ~ 1 /m h rather than power 
of A. We therefore need to extend ([5]i to second order in 
the SCET expansion parameter A. The basis of SCETi 
heavy-light current operators at 0(A 2 ) contains 15 dif¬ 
ferent operator structures, which are listed in ETft . A 
few representative examples are 

Jf = W c ) s iYj'j(x ± D s h v ) 

Jf = ^W c ) Sl iY ± [W ] c iDl c W c \ Sl Tjh v (32) 

and 
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J[ 2 = Qw c ) Sl [[wliD^ c w c ] S2 AwliDl c w c ] S3 }r^ v 


j(2) 

'14 


QW c \Tjh v 


Qw c ) Sl f-r r (wU) s , 


(33) 


where F ) = {1, ys, y„ ) denotes a basis of the four independent Dirac matrices between £ and /z,, . The last two operators 
contain collinear field products at three independent light-cone positions x + Sjti + , i = 1,2,3. The momentum-space 
coefficient functions of these currents therefore depend on two independent variables that describe the distribution of 
the total light-cone momentum among the three collinear fields. These operators describe an effective heavy-to-light 
transition vertex with two additional transverse hard-collinear gluons (/ 12 ), or an additional hard-collinear quark- 
antiquark pair (J 14 ) plus any number of n+ ■ A c gluon fields from the hard-collinear Wilson line W c . The soft fields 
including h v are multipole-expanded and depend only on x p . At this point, we can write the representation of the QCD 
current in SCETi to any order in the A expansion in the form 

OAF, Q)(x) = e- im » vx J] Cf-Cs 1 ,..., s n ) ® Jf\s u ..., s n ; x), (34) 

j.k 


where the ® stands for a convolution over a set of dimen¬ 
sionless variables .?,■ = .v,;;z/,. The superscript Prefers to 
the scaling of the current operator with A relative to the 
leading-power currents. The subscript j enumerates the 
effective currents at a given order in A. The two-point 
correlator ( fTT) of QCD currents is accordingly written 
as 

T» v - H jr (s 1 ,..., s n ) ® tf^h (35) 


quarks read 

4° = K-n-WcgF^W^t, (38) 

4? = q s W)iI/) ±c ^-^ ±c W c q s . (39) 

At second order there are several terms. Examples are 

= (40) 


The hard function /7 ;; . is a product of SCET Wilson co¬ 
efficients C <k) C ( t! and T eS ,'^ v is a correlator of two of the 
i] 1 j jj 

above SCETi currents. 

The effective correlators are computed with the 
SCET] Lagrangian including the power-suppressed in¬ 
teractions up to 0(A 2 ), which have been derived in 
1 H ED. The power-suppressed SCETi Lagrangian 
terms are dealt with in the interaction picture. Writing 

£sce Ti = £ {0) + £ l) + £ {2> + ..., (36) 

the possible time-ordered products that build up the 
hadronic tensor at 0(A 2 ) are 

a) + sym. 

b) 

c ) 7 ( 0 ) 7^ 1 ) -C <1) + sym. 

d) 

e) / 0) / 0) X {1) X (1) . (37) 

The 0{A) suppressed interactions in the effective La¬ 
grangian involving hard-collinear (£) and soft (q s ) 


and 


r( 2) 

■^-HQET.mag 


r 

^mag 

Am h 


Kcr^gF^K, 


(41) 


where C mag A- 1 represents the renormalization of the 
chromomagnetic interaction by hard quantum fluctu¬ 
ations. The complete list of relevant terms in given 

in BP . 


It is worth noting that the interaction in (|39|) de¬ 


scribes the process q c —> g c + q s , that is, the radialion of 
a soft quark from an energetic quark, which is thereby 
seen to be power-suppressed relative to the emission of 
soft gluons, q c —> q c + g s . despite being power- 
suppressed, is a crucial term for hard-spectator scatter¬ 
ing at leading order in the 1 /zzzz, expansion discussed in 
previous sections, since it turns the soft spectator quark 
in the B meson into the energetic quark that forms the 
constituent of an energetic light meson. 4 > on the 
other hand, describes the correction to the soft-gluon 
emission process q c —> q c + g s , and hence the leading 
correction to the usual soft (“eikonal”) approximation. 
Due to the vanishing of odd powers in the A expansion, 
even the second-order corrections in A are needed to 
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obtain the A /nib term to the semi-leptonic decay rate, 
which includes a time-ordered product with two inser¬ 
tions of in (37 i. 

By construction the time-ordered products (37 1 are 
evaluated with the leading-order Lagrangian We 
recall that when the hard-collinear fields are redefined 
according to m 

f = Y^°\ A c = YA^Y', W c = YWf ] Y\ (42) 
where Y is the soft Wilson line of «_ • A s , the soft and 


collinear fields are decoupled in ,£ (0) . The B meson 
by definition contains no collinear degrees of freedom, 
meaning that the B-meson state is represented as the 
tensor product |B) ® |0), where the first (second) factor 
refers to the soft (hard-collinear) Hilbert space. It fol¬ 
lows that the matrix element of any SCETi current cor¬ 
relator 7’ eff /iy , including in general time-ordered prod¬ 
ucts with sub-leading interactions from the Lagrangian, 
can be written in the factored form 


(B|7’ eff (s 1 ,..., S n )\B) = i J d 4 xd 4 y ... e i(mbV ~ q)x 

x (0|[hard - collinear fields]|0)(si 


The additional integrals over cl 4 y ... are related to inser¬ 
tions of the power-suppressed Lagrangian terms. The 
soft matrix element depends only on x+ = n+x,y+ 
n . y,..., so the integrations over transverse positions 
and the n components can be lumped into the defini¬ 
tion of the collinear factor. The soft and collinear ma¬ 
trix elements are then linked by multiple convolutions 

over the light-cone variables x+,y+, _ Having carried 

out these steps, and using momentum-space representa¬ 
tions of the hard coefficients, jet-functions, and shape- 
functions, a generic term in the factorization formula is 
a sum of convolutions 

T = Y i m« u ...«i)*Jiuu . . . , U(\ (J)\, . . . 

(44) 


(B|/z v [soft fields]/!,,|B)(x_,y_,...) 

,...,s„;x,y ,...). (43) 


to any order in the A expansion. The convolution vari¬ 
ables iij = n+ ■ Pi I mb are the fractions of longitudinal 
momentum n+ ■ pi carried by the hard-collinear fields 
in SCETi. The convolution variable oj, is conjugate to 
n+ ■ Xi, and corresponds to n ■ L, where k, is the (outgo¬ 
ing) momentum of a soft field. 

The jet functions are defined as the vacuum matrix 
elements of hard-collinear fields as they appear in ( f43] >. 
Since the hard-collinear virtuality ;;z/,A » A 2 , the jet 
functions can be computed in perturbation theory. On 
the other hand, the soft B-meson matrix elements in ( [43] ) 
define non-perturbative functions. At any given order 
in the A expansion, the soft functions are obtained by 
stripping the hard-collinear fields (since they define the 
jet functions) from the time-ordered product terms. Up 
to 0(A 2 ), one is left with the B-meson matrix elements 


(B\(hvY)(x-)a a (Y 4 h v )(0) b p\B), (B\(h v Y)(x-) aa (Y 4 iiy‘Y)(z-) c d(Y 4 h v )(0)b/3\B), 

(B\(h v Y)(x-) aa (Y 4 iD^ ± Y)(zi-)cd(Y 4 iD v s± Y)(z 2 -) ef (Y 4 h v )(0)b/i\B), 

(B\(h v Y)(x-) aa (Y^h v )(0)bp(q s Y)(zi-)cy(Y^q s )(z2-)ds\B), (45) 


where colour (Latin) and spinor (Greek) indices have been made explicit. These matrix elements are then decomposed 
in colour and spin into scalar functions, which depend on one, two or three to variables. It follows that to order 1 //;(/,, 
but to arbitrary order in a s in the coefficient functions, the semi-inclusive semi-leptonic differential decay rate depends 
on a large number of multi-local shape-functions. An investigation of the structure of the convolution shows that at 
this order the full-QCD hadronic tensor has the factorized structure 


T 


H ■ J(co) ® S (to) 
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+ 2 #(«1 , U 2 ) ®ST(U\,U2'i w) ® 5 (<A) + ^ //(m) ® <Wl, W 2 ) ® 5 (o»l, W 2 ) 

• 6^2, ^ 3 ) ® 5 (6 l>i, 6l>2 , <^> 3 ) + . . . , (46) 











Figure 3: Tree diagrams contributing to the current correlators 7’ eff, ^ v . 
Not shown are diagrams that vanish when n+A c = 0, n-A s = 0, or 
are symmetric to those shown. External lines (except for the dashed 
current insertion) are soft. The double lines refer to the heavy-quark 
field h v . Figure from ED 


where the ellipses denote 1 /m 2 b terms not considered 
here, and for each term the dependence on the con¬ 
volution variables of the most complicated structure is 
shown. The first line represents the leading order in 
the A/mi, expansion. It depends only on a single func¬ 
tion of a single variable. It is evident that the A /nit, 
correction is much more complicated. It is unlikely 
that multi-dimensional soft functions can ever be re¬ 
liably determined from data or be computed by non- 
perturbative methods, so some amount of modelling of 
these functions will be required. The renormalization- 
group equations of the effective theory allow large loga¬ 
rithms In nib /A to be summed to all orders. The anoma¬ 
lous dimensions are themselves multi-dimensional in¬ 
tegration kernels. There is a considerable number of 
them, given the number of different terms, which may 
all mix under renormalization. The problem of sum¬ 
ming logarithms to all orders at <9(1 /mb) has therefore 
so far not been addressed. 

Fortunately, the structure of the result is much sim¬ 
pler than d46l) in the tree approximation as was shown 
in eg eot Figure [3] displays the relevant tree dia¬ 
grams. The reason for this simplification is that tree 
diagrams can only have cuts with a single internal hard- 
collinear particle, and hence there are no convolutions in 
hard-collinear variables Uj. The convolutions of multi¬ 
dimensional soft functions with the tree-level jet func¬ 
tions can be lumped into a set of unknown functions of 


a single variable. An interesting aspect of this analysis 
is that there is a contribution from a four-quark HQET 
operator of the form h Y h v q s q s to the differential semi- 
leptonic decay rate at 0(1/nib) (see the middle diagram 
in the last row of Figure |3j, whereas in the total decay 
rate such contributions are at least of order 1 /m b . For 
a given flavour of q s , these operators contribute differ¬ 
ently to the decay of the charged and the neutral B me¬ 
son, primarily in the region of small n ■ P, where P is 
the total momentum of the hadronic final state X u . For 
a numerical estimate, which requires some modelling, I 
refer to 15T1 . 

Note that the operator T 1 "' in (31 1 appears in similar 
form in many other processes. The SCET formalism 
described here can be adapted to address power correc¬ 
tions in the process e + e —» two jets, in semi-inclusive 
deep-inelastic scattering (DIS) e~ + p —> e~ + X, when 
the hadronic state X is collimated to form a single jet, or 
Drell-Yan production near the threshold. For DIS, the 
analogy is especially clear, since one only needs to re¬ 
place the /(-meson by a proton, and the weak current by 
the electromagnetic current. SCET provides the expres¬ 
sions for the sub-leading interactions and hard vertices 
that are needed to extend the standard factorization the¬ 
orems based on the eikonal approximation and the fac¬ 
toring of jet functions to sub-leading power. 


3. Radiative corrections 

Significant efforts have been dedicated over the past 
ten years to push factorization calculations to 0(a 2 s ) in 
the perturbative expansion. This implies two-loop cal¬ 
culations of the matching coefficient of the SCETi cur¬ 
rent fh Y and the form-factor kernel T l for hadronic two- 
body decays, and one-loop calculations of the hard- and 
hard-collinear matching coefficients in the spectator¬ 
scattering terms. In this section I review these efforts 
and summarize a few results. 

3.1. Form-factor relations 

The matching coefficients relevant to the factorization 
of the heavy-to-light form factors can be read off from 
( [l 6 | . At 0(a 2 ) one needs the one-loop correc¬ 
tions to Cf 1 ’, the coefficient functions of the f c A ± hch r 
operators, and to /, defined by the non-local operator 
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Figure 4: Representative diagrams for the two-loop contribution to C/, 
and the one-loop contribution to cj 51) and J (from left to right). 


matching equation ( TOk. The two pieces were obtained 
in mm, and 0311341 . respectively. In addition, one 
needs the two-loop coefficient functions C, of the g c h v 
operators, which have been calculated in Il52ll53ll54ll55l 
for the vector and axial-vector QCD current and in |56l 
for the remaining (pseudo) scalar and tensor currents. 
Since all quantities involved have operator definitions, 
these are fairly standard loop calculations with special 
vertices derived from the form of the operators and the 
SCET Lagrangians. Figure [4] shows a representative di¬ 
agram for each of the three calculations. 

With these results the symmetry-breaking corrections 
to the heavy-to-light form factor relations Il23l that were 
first computed at 0(a s ) in l24l . can be extended to 
0{cr s ). This has interesting applications to exclusive 
and semi-inclusive semi-leptonic and radiative B decays 
j56[, of which a few will be mentioned here. 

Two examples of form-factor ratios, which simplify 


in the heavy-quark limit, are 


Ki(E) 

fti(E) 


m B f T (E ) 
m B + m P f+(E) ’ 
m B + my T\(E) 
m B V(E) ' 


(47) 


We can write the first in the form 


K\(E) = R t (E) + 


f 


(It C { M\t,E) 


Z p (t,E) 

ME) ’ 


(48) 


where Rp = Cf T /Cf is the ratio of the two-loop match¬ 
ing coefficients for the tensor and vector currents for 
the B —> P (pseudo-scalar) transition, C'f \t,E) = 
Cff(r,E) - C ( j X \t,E) R t (E), and 3 p (t,E) is given in 
( fl5l >. A similar result holds for the second ratio RJE) 
that applies to B —» V (vector) form factors. We recall 
that factorization as discussed here is valid when the en¬ 
ergy of the outgoing light meson satisfies E » A. In this 
limit, the ratios above approach the value 1 plus per¬ 
turbative corrections comprised in R(E ), which are in¬ 
dependent of any non-perturbative parameters, and the 
spectator-scattering correction. Since this term begins 
at 0(a s ), the sensitivity to non-perturbative form fac¬ 
tors and light-cone distribution amplitudes is reduced to 
0(a s ) in ratios. 

To be specific, consider the point of maximal recoil, 
i.e. E = m B /2 or q * 2 3 = 0, at which we find 


f?:r(£max) - 1 + 


47T 




[3 3 J 

l 47T ) 


100 200 A(2 922 

■—EpL v +—L, + 6L;-—L v 


16 10 4 952 , 8047 128 2 , „ 

-—£(3) + — 7T 4 - —7 r + —— + —^ln2 

3 3 27 162 27 


(49) 


with L m = In Qr/mfy, L v = In (v 2 /mj ) ). The scale v is used to distinguish the scale dependence of the QCD tensor 
current from the residual scale dependence of the truncated perturbative expansion. Numerically, the form factor ratio 
7?i(£max) is (setting v = = m h ) |56] 


^l(Emax) = 1 + 


0.046 (NLO) +0.015 (NNLO) 


(El) 


-0.160 |l + 0.524 (NLO spec.) - 0.002 (^ g ) 


= 0.817. 


(50) 


In this expression we separated the symmetry- 
conserving term (first number, equal to 1), the correc¬ 
tion to Rj (remaining terms in the first line) and the 
spectator-scattering contribution (second line, including 
a small correction dj og from renormalization-group sum¬ 


mation |[34l ), and within these the 0(a 2 s ) contributions 
“NNLO” and “NLOspec.”. 

We observe that these 0{a 2 s ) corrections are 30% 
(first line of ( |50] >) and 50% (second line) of the 0(a s ) 
ones, and result in an overall reduction of 'R i (£,’ rn ax ) 
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Figure 5: The matching coefficient cj(w, fi) / c 9 fu, /i) as a function of 
u (related to the di-lepton invariant mass q 2 = (1 - u)m 2 ) in the one- 
loop (dashed) and two-loop (solid) approximation. The blue/dark grey 
curves refer to fi = mt, =4.8 GeV, and the orange/light grey ones to 
/i = 1.5 GeV. Figure from EQ 


by 20% relative to the symmetry limit. The Rj and 
spectator-scattering corrections have opposite sign, but 
the latter are larger and determine the sign of the devi¬ 
ation from the symmetry limit. The same observations 
hold for ?? 2 (£maxX but in this case the sign of the two 
terms is opposite and one finds 7?2(£max) = 1.067. The 
normalization of the spectator-scattering term (the num¬ 
ber -0.160 in ([50])) involves the parameter combination 


9 fJ B 

m b f+ 7r (0)A B ’ 


(51) 


which captures most of the dependence on non- 
perturbative parameters of spectator scattering and the 
form-factor ratio as a whole, and supplies the by far 
dominant source of quantifiable theoretical uncertainty. 
In addition there are 0(h/mb) corrections to the form- 
factor factorization formula (J9]l. 

There are also QCD sum-rule calculations of the 
QCD form factors, which have different theoretical 
systematics, and implicitly include some 0(A/nib ) 
power corrections. These calculations give mm 
R\(E mm ) = 0.955 and 7? 2 (£max) = 0.947. For the 
tensor-to-vector ratio R 2 , one notices that the sign of 
the symmetry-breaking correction is opposite for these 
two methods. The origin of this discrepancy is not un¬ 
derstood. Before drawing conclusions on the sum-rule 
method or the size of power corrections, however, a 
dedicated analysis of form-factor ratios (rather than the 
form factors themselves) with correlated theoretical un¬ 
certainties should be performed within the QCD sum- 
rule method. 


3.2. Semi-inclusive B —> X s ( + { decay 

The two-loop correction to the QCD current match¬ 
ing to SCETi has an interesting application to semi- 
inclusive B —* X s { + {~ decay. As discussed in Sec¬ 
tion [23] the semi-inclusive rate factorizes in the form 

dT [0] = h [0] xJ®S , (52) 


at leading order in the 1 /nib expansion with hard func¬ 
tions / 2 101 that can be expressed in terms of the coeffi¬ 
cients C; Of £.r,/h- 

Here I focus on the differential (in q 2 , the ( + t~ invari¬ 
ant mass) forward-backward (FB) asymmetry in the an¬ 
gle between the positively charged lepton and the B me¬ 
son in the centre-of-mass frame of the di-lepton pair in 
the presence of an upper limit on the invariant mass nix 
of the hadronic final state I 59I160B . It is well-known that 
due to the interplay of electroweak penguin operators 
mediating b —> s{ + £~ and the electromagnetic dipole 
operator mediating b —> s~y(—*£ + £~), the FB asymmetry 
exhibits a zero. In the presence of an invariant-mass cut 
the location q^ of the zero is given approximately 
by I56]| 


% 


_ Re[C“ cl (^)] c](hq) 

2 m b (niB - (p*)) “ Re [C“ cl (^)] c > 0 ) 


(53) 


with uq =1 ql/(m h (m B - (p +») and </?+> w p+ cut /2 . 
This result depends on ;«™‘ through the corresponding 
cut on the light-cone momentum component /?£ cut . 

The first factor on the right-hand side of (53 1 is pri¬ 
marily related to the Wilson coefficients of the elec¬ 
troweak penguin and electromagnetic dipole operators 
in the weak effective Lagrangian. The second factor 
Cj(m)/Cj(m) provides a modification from the matching 
to SCETi and can be expressed in terms of the coeffi¬ 
cients C, known to the two-loop order. The dependence 
of this factor on the di-lepton invariant mass is shown at 
LO, NLO and NNLO in Figure[5]for two different renor¬ 
malization scales. The impact of the NLO correction is 
to shift q^ by -2.2%. The size of the NNLO correction 
is significant, in fact larger than the NLO correction, and 
shifts the zero by another -3%. Including an estimate 
of power-suppressed effects. Ref. ff56l finds 


< 7 0 2 = [(3.34... 3.40)!°;g]GeV 2 (54) 


for m = (2.0... 1.8) GeV. This value of the asymme¬ 
try zero in semi-inclusive b —> s( + t~ decay is signifi¬ 
cantly smaller than for the exclusive case (see below), 
where spectator scattering is responsible for a positive 
shift as is the fact that in this case ( p J) = 0 in (53 1 . 
On the other hand the semi-inclusive zero is in the same 
region as in the inclusive case El. 
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Figure 6: Representative diagrams for the one-loop contribution to 
H n (left), and the two-loop contribution to T l (right). 


3.3. Tree-dominated charmless hadronic decays 


The matching coefficients relevant to the factoriza¬ 
tion of the charmless, hadronic, two-body decays can 
be read off from ( pT) . At 0(a 2 s ) one needs the 

one-loop correction to H n , the coefficient function of 
the (xx)(£cA±hchv) spectator-scattering operator C? 11 , and 
the two-loop correction to T l related to the (xx)(£,h v ) 
operator O 1 . The hard-collinear function J is the same 
that enters the form factor. The computation of H 11 at 
one-loop was performed in Il38ll for the so-called tree- 
operators Q\ 2 in the weak effective Lagrangian. Two 
independent calculations can be found in lf62l[63 |. The 
two-loop matching of £> 1,2 to the SCETi four-quark op¬ 
erator, which defines T l , was first done for the imagi¬ 
nary part of the amplitude l64l and then in two indepen¬ 
dent calculations l65l[66 l for the full amplitude, finding 
complete agreement. These results together extend the 
0(a s ) result of mini to the next order for pure tree- 
operator decays such as B —> n n°. Figure [6] shows a 
representative diagram for each of the two calculations. 

The calculations are more involved than those for the 
form factors, since the presence of the second light me¬ 
son adds another variable to the functional dependence, 
such that H u depends on two momentum fractions, see 


(201. Furthermore, the four-fermion operators consist of 


O n = y7^7± ■ • -7± ‘C 1 - 7s) ® yO - 


with n = 2,3,4. All these operators are evanescent, 
i.e. vanish in four dimensions. They disappear from 
the final result, since we shall renormalize them such 
that their infrared-finite matrix elements vanish, but they 
must be kept in intermediate steps, hence the matching 
equation ( fT7] > has to be extended to include all four op¬ 
erators on the right-hand side. 

Evanescent operators appear already at tree level. In 
this approximation the right-insertion matrix element of 


two strings of fermion lines between which gluons can 
be exchanged. When factorization is done in dimen¬ 
sional regularization, one has to deal with the renor¬ 
malization of operators that vanish in four dimensions 
(“evanescent operators”), whose matrix elements must 
be renormalized to zero. Related to this is the issue of 
consistent Fierz transformations in d dimensions. The 
transition b —> u(itd), which is effected by the local op¬ 
erators Q\ 2 and where the bracket indicates the fermion 
fields which are contracted, can contribute to the n'n 
final state in the form ( q s u)(iid ) (called “right insertion”) 
as well as to the 7r°7r° final state in the form ( q s d)(uu ) 
(called “wrong insertion”), where q s = d is the spectator 
quark in the B c i meson. In the second case a Fierz rear¬ 
rangement of the original fermion strings is required to 
match to the SCETi operators ( p~8j ), | [T9) . The simplest 
way to deal with this difficulty is to add the difference 
of Fierz-equivalent operators to the list of evanescent 
operators to be renormalized to zero Il38l (66 1. In the 
following I briefly elaborate on the issue of evanescent 
operators and the subtractions required to arrive at the 
renormalized coefficient functions. 

In case of H u (u , v), one has to calculate the one-loop 
five-point amplitudes b —> q C 2 <ic 2 <lc\gc\ associated with 
the matrix elements (q(qi)q(q 2 )q(p' l )g(p' 2 )\Qi. 2 \Kp)), 
where the collinear light particles carry momentum q\ = 
umbii+l2 , q 2 — umbn+12 in collinear-2 direction, and 
p\ = vmbn-l 2, p\ - vmbn-12 in the other. When one 
strips the SCETj operator 0 ll (t, s) off all its fields, and 
represents it only by its Dirac stmcture, 

Oi ee yd- 75 ) ® y(1-75)7^, (55) 

one finds that one-loop right-insertion amplitude also 
contains the operators 


75)7±a, ■ ■ ■ 7 - 1071-1 


(56) 


Q 2 is given by 


(22)i it - — \ -{ 0 \) - -{O2) 

JSI c \u uu 


(57) 


The subscript “nf” (for “non-factorizable”) means that 
the “factorizable” diagrams with no lines connecting 
the collinear-2 and collinear-1 sectors are omitted, since 
they belong to the T ] O l term in (17 1 by virtue of the 


definition (fl8|. While one can simply set (O 2 ) = 0 in 
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the above equation, since no 1 /e poles are present at 
tree level, the appearance of an evanescent operator in 
the ^-dimensional tree level matrix element implies that 
one must compute the mixing of Os into 0\ in the 1 - 
loop calculation. The renormalized coefficient function 
of the physical (non-evanescent) operator O 11 can then 
be written as 


H m) = A' 1 ' - A' 0) * Z ( H - Af * Z™ 


(58) 


where A® denotes the coefficient of O n of the non- 
factorizable contribution to the ultraviolet-renormalized 
QCD /f-loop b —> qciqaqcigcx amplitude. The first 
term is the one-loop amplitude, which is infrared di¬ 
vergent. The next term is the SCETi MS counterterm 
Zjj for the physical operator times its tree coefficient 
A®, which subtracts this divergence. The third term 
arises from the mixing of the evanescent operator O 2 
into 0\ at one-loop. The finite renormalization con¬ 
stant Z \= -M**’ 0 * can be computed from the mixing 
matrix element with an infrared regularization different 
from the dimensional one (for instance, off-shell). Fi¬ 
nally, the last term arises from the fact that we defined 
( p~ 8 1 to reproduce the full QCD form factor. The stars 
in ( 58]) remind us that the amplitudes depend on mo¬ 
mentum fractions, so the renormalization “constants” 
are actually convolution kernels. For the wrong inser¬ 
tion coefficient function, the expression is slightly more 
complicated and contains an additional term 



It is finite and independent of the infrared regulator and 
ensures that the difference of Fierz-equivalent SCETi 
operators is correctly renormalized to zero. Further¬ 
more, due to the Fierz rearrangement there could be a 
contribution 


A«-A« + A^*M<r 

+ A^*« off -M^ ff ), (60) 

from factorizable diagrams, which is not contained in 
the definition of the QCD form factor. However, this 
term vanishes here. 

The calculation of T ] (u) amounts to evaluating the 
two-loop on-shell Q, matrix element of the transition 
b(p) —» q\(p')q 2 (uq)q 3 (uq) with kinematics p — p' + q, 
p- = m 2 b , p' 2 = q 2 - 0. The ultraviolet and infrared 
finite short-distance coefficient is obtained from 


y 'tO j 
± i 
y(l) 


i(0) 


A ( ! )nf +Z ( 1 ) A ( °\ 


y( 2) 


A (2)nf 

+ 

N 

+ 

Z (2) A ( ? 

+ Z®A ( ! )nf 


ij Ji 

ij j 1 

a 1 1 

+ (—i) 8m (V> A'[ 1)nf 



_ y( 1) 

x i 

+ 

— 

C fc, 

0 

-Z (1) l 

Zj ext\ 


-Z 



(61) 


b> 1 


where now A® denotes the bare A-loop QCD amplitude, 
Zy, the QCD renormalization matrices including 
those for the operators Q , and the bottom mass coun¬ 
terterm, Y a t, the SCETi renormalization factors, and 
1 the one-loop coefficients of evanescent SCETi four- 
quark operators. Again convolutions are implied, for 
instance, if® TfJ' must be interpreted as the convolu¬ 
tion product j' du' H^\u') Y^{u', u). The correspond¬ 
ing equation for the two-loop wrong insertion matrix 
element is significantly more complicated, as the ad¬ 
ditional terms from Fierz rearrangement do not vanish 
any more, and the difference of two-loop SCETi coun¬ 
terterms must be computed ( 66 l . 

The computation of the bare two-loop QCD matrix 
elements A (2,n( is technically the most challenging part 
of ( |6T| l and has been done by reducing tensor integrals 
to scalar integrals by Passarino-Veltman reduction ll67| , 
the reduction of scalar integrals to master integrals us¬ 
ing the Laporta algorithm ll 68 l based on integration-by- 
parts (IBP) identities 169/70], and finally by evaluating 
the master integrals directly, or through Mellin-Barnes 
representations in the more complicated cases. It is con¬ 
venient to compute the convolutions 

6w(l - h)C< 3/2) (2k - 1) Tf\u) (62) 

of the kernels with the first few terms in the Gegenbauer 
expansion 



<Pm(u) = 6 uu 


l+2afCi 3/2) (2«-l) 


(63) 


of the light-meson LCDA to express the final result 
in terms of the first two Gegenbauer moments a^ 2 . 
In ll 66 l fully analytic expressions after integration over 
the Gegenbauer expansion, including the exact depen¬ 
dence on the charm-quark mass, have been obtained. 

Putting these results together, we can proceed to the 
investigation of the so-called topological tree ampli¬ 
tudes with <9(o; 2 ) accuracy. For definiteness, only the 
three B —» nn decay modes will be considered. The 
decay amplitudes are given by 


V2 )71b ji V] — A nn A u 


a 1 + a 2 
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— >7T + 7T~ ~ A-nn | 

a 1 + a\ 

+ A c <3(4 j 

— ■^ljjo_ > 7 r o 7r o — A nn | A u 

a 2 - a\ 

— A c Q '4 j 


where A nn = iG F m\f„f+”(0)l V2 and A p = V pb V* d . 
The dominant amplitudes governing the set of tut 
amplitudes are the so-called colour-allowed (colour- 
suppressed) tree amplitude a\ (nn) (a^Jnn)) and the non¬ 
singlet QCD penguin amplitudes a P A (nn). The equa¬ 


tion does not show some smaller amplitudes that are 
taken into account in the numerical evaluation of the 
branching fractions below. (Full expressions for the 
amplitudes are given in IfTTl .) The normalization of 
the topological tree amplitudes is an = 1 and a 2 = 
1/3 in the approximation by tree-level Vk-hoson ex¬ 
change, when renormalization-group evolution of the 
Qi from the electroweak scale to m b is neglected. The 
NNLO QCD factorization results for the tui final states 
read ll66l 


a 1 (tui) 


a 2 (n-7r) 


1.009 + [0.023 + 0.010 /] NL0 + [0.026 + 0.028 /]nnlo 


• sp 


{[0.014] LOsp + [0.034 + 0.027/] NLOsp + [0.008] tW 3 J 


0.445 

0.220 - [0.179 + 0.077/Jnlq - [0.031 + 0.050/] NNL0 


' sp 


0.445 J 


{[0.114] LOsp + [0.049 + 0.051/Wp + [0.067] tw3 j 


1-000-U069 + (0-01 !-o:o!o)d (65) 

0.240^ + (-0.077!°;“ 5 8 )/. ( 66 ) 


The first line of each of the two equations corresponds 
to the contribution from the operators ixx)QcK) with 
coefficient function T 1 , while the second line accounts 
for spectator scattering at 0(a s ) (“LOsp”) and 0((y 2 ) 
(“NLOsp”) and a certain power correction (“tw3”). The 
numerical result sums all contributions and provides an 
estimate of the theoretical uncertainties from the param¬ 
eter variations as detailed in lt 66 l . 


The amplitudes aio are scale- and scheme- 
independent physical quantities. Figure [7] shows the 
residual renormalization scale dependence at LO, NLO 
and NNLO of the vertex correction T l to the colour- 
suppressed amplitude a 2 . The NLO contributions to 
this amplitude are large, because they are proportional 
to the tree operator with large Wilson coefficient in the 
weak effective Lagrangian. The scale dependence sta¬ 
bilizes for the real part at NNLO, while the absorptive 
part receives another sizable correction, approximately 
25% of the leading-order real part, and hence its scale- 
dependence is barely reduced. By comparing the theo¬ 
retical uncertainty of the full expressions in (65 1 , (661 
above to the one from the vertex correction alone, see 
Figure [7] we see that it arises primarily from spectator 
scattering. The main contributors to the uncertainty are 
the parameter combination 0 - which appears as an 
overall normalization factor of the spectator-scattering 
term, the second Gegenbauer moment a/(2 GeV), and 
the “tw3” estimate the power correction. 


We observe that the vertex and spectator-scattering 
corrections come with opposite sign, both at NLO and 
NNLO, and separately for the real and imaginary part. 
This seems to be a general pattern that was also ob¬ 
served for the form factor ratios discussed in the previ¬ 
ous subsection. Here it is somewhat unfortunate, since 
an enhancement rather than a cancellation in the colour- 
suppressed tree amplitude would have been welcome in 
view of the trend indicated by experimental data, as will 
be seen below. 


The colour-suppressed tree amplitude a 2 {mi) in par¬ 
ticular exhibits an interesting structure. It starts out with 
a positive real value 0.220 at LO. After adding the per¬ 
turbative corrections to the four-quark vertex (the first 
line of ( 661 ), it is found to be almost purely imaginary, 
0.01 - 0.13/. Then the spectator-scattering mechanism 
regenerates a real part of roughly the original size and 
cancels part of the strong phase. The net result of this 
is that the colour-suppressed tree amplitude can become 
sizable in QCD factorization when r sp is large, but since 
this enhances the cancellation of the imaginary part, 
one cannot have both, a large magnitude and a large 
strong phase. In essence, the dynamics of the colour- 
suppressed tree amplitude is completely governed by 
quantum effects, and the theoretical uncertainty is cor¬ 
respondingly large. In comparison, the colour-allowed 
tree amplitude ai(7m) is rather stable against radiative 
corrections, and never deviates by a large amount from 
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Figure 7: Dependence of the topological tree amplitude a on 
the hard scale \i (vertex correction T l only). The dotted, dashed and 
solid lines refer to the theoretical predictions at LO, NLO and NNLO, 
respectively. Figure from ED- 


its LO estimate. 

Strictly speaking, only the charged B-meson tree- 
operator decays to pions and longitudinally polarized 
rho mesons can presently be predicted with complete 
NNLO accuracy, since the penguin amplitude a p that 
appears in the other two modes of <(64j is not yet com¬ 
pletely known at 0(a 2 s ) (see next subsection). By us¬ 
ing the best available approximations for a p the full 
set of 7 T 7 T, npi and PlPl final states has been ana¬ 
lyzed in ( 66 , 721. Table [2] shows an excerpt of the 
predicted CP-averaged branching fractions compared 
to present measurements compiled from BABAR and 
BELLE (including CDF and LHCb for 7t + n j data as 
given in urn Within uncertainties the agreement be¬ 
tween theory and measurement is very good for the fi¬ 
nal states with a substantial contribution from a i. The 
colour-suppressed amplitude ai appears to be predicted 
too small, although the comparison recently improved 
due to a (preliminary) revision of the previous BELLE 
measurement. The uncertainties in Tableware strongly 
correlated among the three decay modes. Further infor¬ 
mation can be brought to light by considering ratios of 
branching fractions sensitive to certain parameters and 
amplitudes. I refer to (6611721 for the details of this anal- 


Theory Experiment 

B -> 5 . 82 “;« (*) 5 . 48 *$ 

B° d —> tTtC (*) 5.10 ±0.19 

B° n°n° 0.63^:1H:« 1-17 ± 013 

Table 2: CP-averaged branching fractions (BrAv) in units of 10 ( ' of 
B —» nn decays. The first error on a quantity comes from the CKM pa¬ 
rameters, while the second one stems from all other parameters added 
in quadrature. Theory corresponds to Theory II in (66) and adopts 
the values 0) = 0.23 ± 0.03, T B (lGeV) = (0.20*>-°j’) GeV for 
the pion form factor and B-meson LCDA parameter, respectively. In 
order to focus on the hadronic uncertainty, the ranges given refer to 
the uncertainty of BrAv( B —> f)/\V u b\ 2 - Ranges marked by an aster¬ 
isk (★) handle the dependence on the heavy-to-light form factors in a 
similar way and refer to BrAv(B —» f)/(\V u b\ 2 f? n (0) 2 ). 


ysis and results for final states containing rho mesons. 

3.4. Penguin-dominated charmless hadronic decays 

Direct CP asymmetries require the interference of 
amplitudes with different CKM factors, A u and A, in 
as well as a phase difference of the corresponding 
hadronic matrix elements. This necessitates the compu¬ 
tation of the matrix elements of the so-called penguin 
operators Q 3-10 and dipole operators Qn y , Qg g in the 
weak effective Lagrangian. 

At 0(ar s ) the one-loop coefficients Hf and two-loop 
coefficients Tj for < 23-10 (one loop less for Qi y and Qg g ) 
are needed. They involve, in addition to diagram topolo¬ 
gies identical to those relevant to Q\ 2 such as shown in 
Figure[ 6 ] the “penguin contractions” of same-flavour qq 
fields from the operators Q,. On the one hand, the pen¬ 
guin contractions are topologically simpler than the ver¬ 
tex diagrams, on the other hand the presence of an inter¬ 
nal massive charm or bottom quark loop introduces an¬ 
other dimensionless ratio of scales, which complicates 
the analytic calculation of the master integrals. The cal¬ 
culation of the hard spectator-scattering kernels for the 
non-singlet QCD and the colour-allowed and colour- 
suppressed electroweak penguin amplitudes has been 
performed in l74l . but the two-loop calculation of T 1 has 
not yet been completed. Technical results on the rele¬ 
vant master integrals have already appeared (75). Thus, 
at present, the 0{a s ) results of J5J [IT] [TTl for the pen¬ 
guin amplitudes and direct CP violation cannot be ex¬ 
tended to the next order. 

I therefore refrain from a detailed discussion of the 
numerical impact of the spectator-scattering contribu¬ 
tion (see l74l ) except for one remark. In principle, the 
one-loop correction to a r 4 can be expected to be rather 
large, since there is a contribution from the tree operator 
Q 1 with coefficient Ci ~ 1 ten times larger than the co- 
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efficients C 3_6 of the QCD penguin operators. However, 
there is an almost complete, possibly accidental cancel¬ 
lation between the contributions from Q\ with colour 
factors Cp — 4/3 and Ca = 3, respectively. 

4. Phenomenology 

Since the focus of this article is on the development 
of factorization and radiative corrections, and since the 
phenomenology of large sets of final states is not yet 
possible at NNLO in the absence of the calculation of 
at least the QCD penguin amplitude with this accuracy, 
this section will be brief, drawing in part on a sum¬ 
mary prepared for j76l . highlighting what in my judge¬ 
ment are interesting conclusions from and applications 
of 0(a s ) (NLO) calculations. 

4.1. Charmless, hadronic two-body decays 

Charmless, hadronic two-body decays provide by far 
the most observables due to the sheer number of pos¬ 
sible final states. Two-body here includes final states 
with unstable particles like kaons, rho mesons etc. as 
long as they can be clearly identified as sharp reso¬ 
nances. NLO QCD factorization results are available for 
a variety of complete sets of final states: two pseudo¬ 
scalar mesons (PP) and pseudo-scalar plus one vector 
meson (PV) 11711 . two vector mesons (VV) 1771 . final 
states with a scalar meson ED, an axial-vector meson 
1791®, and a tensor meson liffl . 

A key issue for phenomenology is the treatment of 
power corrections, since factorization as embodied by 
o is not expected to hold at sub-leading order in 
l/mt,. Some power corrections related to scalar cur¬ 
rents, which appear after Fierz rearrangements, are en¬ 
hanced by large (“chirally enhanced”) factors such as 
m^/((m u + nid) A). Some corrections of this type, in par¬ 
ticular those related to scalar penguin amplitudes never¬ 
theless appear to be calculable and turn out to be impor¬ 
tant numerically. On the other hand, attempts to com¬ 
pute sub-leading power corrections to hard spectator¬ 
scattering in perturbation theory usually result in in¬ 
frared divergences, which signal the breakdown of fac¬ 
torization. These effects are then estimated and included 
into the error budget, see the contributions marked 
“tw3” in ( [65] ), ( |66] >. All weak annihilation contribu¬ 
tions belong to this class of effects and often constitute 
the dominant source of theoretical error, in particular 
for the direct CP asymmetries. Factorization as above 
applies to pseudo-scalar flavour-non-singlet final states 
and to the longitudinal polarization amplitudes for vec¬ 
tor mesons. Final states with rj , //' require additional 


considerations, but can be included f82l . The trans¬ 
verse helicity amplitudes for vector mesons are formally 
power-suppressed but can be sizable ll83ll . and they do 
not factorize in a simple form IfTTl . The description 
of polarization is therefore more model-dependent than 
branching fractions and CP asymmetries. Besides these 
conceptual uncertainties, the lack of precise knowledge 
of quantities such as \V u b\, heavy-to-light form factors, 
and the fi-meson LCDA parameter A B cause a signifi¬ 
cant theoretical uncertainty. 

Much of the phenomenology of the most widely 
considered PP, PV and VV final states is determined 
by the non-singlet QCD penguin amplitude a%. This 
amplitude is certainly underestimated at NLO in a s 
and leading order in the heavy-quark expansion. The 
power-suppressed but chirally-enhanced scalar penguin 
amplitude, and probably a (difficult to disentangle) 
weak annihilation contribution is required to explain the 
penguin-dominated PP final states. While the scalar 
penguin amplitude is calculable, some uncertainty re¬ 
mains. An important observation is the smaller size of 
the PV, VP and VV penguin amplitudes as compared 
to PP final states, which can be inferred from the mea¬ 
sured branching fractions of hadronic b —> s transitions. 
This is a clear indication of the relevance of factoriza¬ 
tion, which predicts this pattern as a consequence of 
the quantum numbers of the operators Qj. If the pen¬ 
guin amplitude were entirely non-perturbative, no pat¬ 
tern of this form would be expected. A similar state¬ 
ment applies to the final states, where factoriza¬ 

tion explains naturally the strikingly large differences 
in branching fractions, including the large rj’K branch¬ 
ing fraction, through the interference of penguin am¬ 
plitudes, although sizable theoretical uncertainties re¬ 
main 182] , The flavour-singlet penguin amplitude seems 
to play a sub-ordinate role in these decays. 

The situation is much less clear for the strong phases 
and direct CP asymmetries. A generic qualitative pre¬ 
diction is that the strong phases are small, since they 
arise through either loop effects ( a s (nib )) or power cor¬ 
rections (A /nib). Enhancements may arise, when the 
leading-order term is suppressed, for instance by small 
Wilson coefficients. This pattern is indeed observed. 
Quantitative predictions have met only partial success. 
The observed direct CP asymmetry in the tt' t: and the 
asymmetry difference in the n°K + , n K 1 final states are 
prominently larger than predicted. A comparison of all 
CP asymmetry results shows a presently ununderstood 
pattern of quantitative agreements and disagreements. 
Since a s (rnb)/n and A/nib are roughly of the same or¬ 
der, it is quite possible that power corrections are C?(l) 
effects relative to the perturbative calculation, prevent- 


M. Beneke, Soft-collinear factorization in B decays (2015) I | 2 /| 


22 


ing a reliable quantitative estimate. However, the di¬ 
rect CP asymmetry calculations are still LO calcula¬ 
tions, contrary to the branching fractions, so the final 
verdict must await the completion of the NLO asymme¬ 
try calculation, which requires the two-loop computa¬ 
tion of a p . Contrary to direct CP asymmetries, the S 
parameter that appears in time-dependent CP asymme¬ 
tries is predicted more reliably, since it does not require 
the computation of a strong phase. This is exploited 
successfully in NLO computations of the difference be¬ 
tween sin 2 \p from b —* s penguin-dominated and from 
b —» ccs tree decays ||84| [85], and the direct determina¬ 
tion of the CKM angle y (a) from time-dependent CP 
asymmetry measurements for the ;r + 7 r~, tF p" and PlPl 
final states ( IfTTl and Figure 1 of ( 86 ]). 

Polarization in B —» VV decays was expected to 
be predominantly longitudinal, since the transverse he- 
licity amplitudes are A/nib suppressed due to the V-A 
structure of the weak interaction and helicity conserva¬ 
tion in short-distance QCD. While this is parametrically 
true (with one exception, see below), closer inspection 
shows that the parametric suppression is hardly real¬ 
ized in practice for the penguin amplitudes (831 . This 
leads to the qualitative prediction (or rather, in this case, 
postdiction) that the longitudinal polarization fraction 
should be close to 1 in tree-dominated decays, but can 
be much less, even less than 0.5, in penguin-dominated 
decays, as is indeed observed. However, quantita¬ 
tive predictions of polarization fractions for penguin- 
dominated decays must be taken with a grain of salt, 
since they rely on model-dependent or universality- 
inspired assumptions of the non-factorizing transverse 
helicity amplitudes (77l . 

There is an exception to the power counting for the 
helicity amplitudes, which arises from a contribution 
of the electromagnetic dipole operator Qi y to the trans¬ 
verse electroweak penguin amplitude o^’p W (87], which 
is enhanced by a factor nib/A relative to the longitudi¬ 
nal amplitude (but proportional to the small electromag¬ 
netic coupling). This effect manifests itself in branch¬ 
ing and polarization fractions of the pK* final states and 
should be seen by comparing precise measurements to 
theoretical predictions and measurements of the related 
PP, PV modes ::K, pK , nK*. In particular, polarization 
measurements would then result in an indirect measure¬ 
ment of the (virtual) photon’s polarization in the elec¬ 
tromagnetic dipole transition, and hence be sensitive to 
the possible presence of a dipole operator with opposite 
chirality of the light quark due to non-standard physics. 



Figure 8: Leading-order diagrams for the B —» yW*(—> tv ) transi¬ 
tion. The left graph shows the leading-power contribution from pho¬ 
ton emission from the up antiquark. Emission from the heavy 6-quark 
(right) is power-suppressed. Figure from (HD. 

4.2. B —> ytv 

The decay B ytv is accessible to factorization 
methods when the energy E y of the photon is large com¬ 
pared to A. Even though the final state does not con¬ 
tain a hadron, the coupling of the on-shell photon to 
soft and collinear quarks is non-perturbative and leads 
to hadronic (0 |... | B) matrix elements of non-local op¬ 
erators. The two tree diagrams relevant to the process 
are shown in Figure [ 8 ] Photon emission from the heavy 
quark (right diagram) is power-suppressed relative to 
the emission from the light quark in the left diagram. 
The differential branching fraction can be expressed in 
complete generality in terms of two independent B —* y 
form factors F Vt A(E y ). 

Soft-collinear factorization properties of the B —» y tv 
amplitude have first been discussed in (891 , and in (90l 
El with SCET methods. At leading order in the l/nib 
expansion, the two form factors coincide and satisfy a 
factorization formula 

F — C ■ J ® <p B . (67) 

In the tree diagram in the left Figure[ 8 ]the hard function 
C is the weak-decay vertex, while the jet function con¬ 
tains the hard-collinear light-quark propagator that con¬ 
nects the VV’-boson and photon vertices. As in other ap¬ 
plications of leading-power factorization, the B-me son 
LCDA enters only through Ab and the logarithmic mo¬ 
ments ( 22 l >. 

The radiative-leptonic decay is the simplest exclusive 
decay, and the /Lineson LCDA is the only hadronic pa¬ 
rameter relevant to this process at leading order in the 
heavy-quark expansion. It therefore provides a unique 
way to determine A s from (future) data. The sensitivity 
to Ab can be seen from the leading-power expression 

FvAEy) = R(E y ,p) ( 68 ) 
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for the B — > y form factors. The radiative correction 
factor R(E 7 ,/j) and the leading power-suppressed ef¬ 
fects have been investigated in (88) . Figure [9] shows the 
branching fraction with a lower cut on the photon en¬ 
ergy as a function of the parameter Ab and confirms the 
strong sensitivity. It can be shown that at sub-leading 
order in A/nib there appears a single unknown “soft” 
form factor similar to the soft form factor in the B —> n 
transition. However, the difference between F v and F A 
is independent of this form factor and can be given en¬ 
tirely in terms of the /I-meson decay constant /#, at least 
at tree level (88). The power-suppressed soft form factor 
limits the accuracy of the Ab determination. A numeri¬ 
cal estimate has been obtained with the QCD sum-rule 


method l92l . 

4.3. Radiative and electroweak penguin decays 

I cannot conclude this overview without a brief men¬ 
tioning of the exclusive radiative and electroweak pen¬ 
guin decays B —> K*y and B —> K (t) t + t~, respectively. 
These are primarily driven by the two loop-induced 
( £€){sb ) operators in the weak effective Lagrangian with 
a non-negligible contribution from the electromagnetic 
dipole operator and four-quark operators. The hadronic 
matrix elements of the former two contributions can be 
parametrized by heavy-to-light form factors, but those 
of the four-quark operators are non-local. Schematically 
the transition matrix element contains two terms. 


{K*tt\H c „\B) = Y, a,(C 


'(') 

7y,9.10 




q 


f- 


d 4 x j** (K*\T(/ em (x)H^ A mB) 


(69) 



Figure 9: Dependence of the partial branching fractions Hr ( li —* 
ytv,E y > E cut) for E ait = 1 GeV (upper band) and 1.7 GeV (lower 
band) on Ab- Figure from (88). 

where the first comprises the electroweak penguin and 
electromagnetic dipole transitions, can be expressed in 
terms of form factors, and is sensitive to other virtual 
particles in the loops that generate these transitions. The 
second is non-local, dominated by QCD effects, includ¬ 
ing charmonium resonances in the b —> s(cc —» It) 
transition. 

Soft-collinear factorization is relevant in two ways in 
the large-recoil region where the kaon energy E K » 
A. First, in the first term the QCD form factors can be 
expressed in terms of only two “soft” form factors 
through the factorization theorem 0 for form factors. 
Second, and more importantly, the complicated second 


term can be calculated and expressed in terms of exactly 
the same form factors, such that the entire amplitude 
( |69| l is expressed as 

= c,f + cPb ® Ti 9 <f>K * + < 5 | — 1 (70) 
\m h j 

for every operator Q, in /7 e fi as shown in [93] [94]| lor the 
radiative decays B —> Vy and in (93) for B —> Vt + A. 
These papers also performed the 0(a s ) calculations. 
Note that NNLO computations for these decay modes 
are significantly harder than for charmless decays, since 
at this order three-loop diagrams must be considered 
for the vertex kernels, and two-loop diagrams for hard- 
spectator scattering. Only the simpler contributions 
from the dipole operators have been calculated up to 
now (95l . The results of (93l [94) were soon generalized 
to isospin asymmetries and the B —> py, ptl 

decays (98) . 

As an example of the uses of factorization I show 
a plot from 193) for the differential forward-backward 
asymmetry already discussed above in the context of 
semi-inclusive hadronic final states. Figure [TO] shows 
the reduction of the theoretical uncertainty at the loca¬ 
tion of the asymmetry zero due to a reduced sensitiv¬ 
ity to the form factors (99l 1 10 01 and a sizable shift of qi 
from q 2 0 = 3.4+J® GeV 2 at LO to q 2 Q = 4.39+°]® GeV 2 at 
NLO. Including an estimate of power corrections to the 
form factors, q = (4.2 ± 0.6) GeV 2 (93l remains a con¬ 
servative estimate of the asymmetry zero location. A 
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Figure 10: Forward-backward asymmetry dA^(B~ —> 

K*~£ + 1~)/dq 2 at next-to-leading order (solid cent line) and leading 
order (dashed) as function of lepton invariant mass q 2 in GeV. The 
band reflects all theoretical uncertainties from parameters and scale 
dependence combined. Figure from ED 

subsequent analysis with updated parameters fl98l gives 

ql[K*°] = 4.36^q3j GeV 2 (71) 

ql[K* + ] = 4.15 ± 0.27 GeV 2 (72) 

for the neutral and charged /I-meson decay separately, 
which replaces the value q 5 = 4.39^j ^ GeV 2 above. 

With the start-up of the LHCb experiment at the hori¬ 
zon and its large anticipated sample of B —» t ] de¬ 
cays, the fully differential angular distribution includ¬ 
ing K* —> Kn became the subject of investigations 
mmnm From the angular coefficients one can con¬ 
struct observables with reduced sensitivity to B-meson 
form factors for all values of q 2 rather than a single 
value as is the case for the forward-backward asym¬ 
metry. This can be done by selecting observables that 
are proportional to only one of the products |£||(c/ 2 )| 2 , 
l^j-P? 2 )! 2 . g±(q 2 )g\\(q 2 ) of soft form factors, and then tak¬ 
ing an appropriate ratio, such that the dependence on £ 
cancels out at leading order in a s . The main theoret¬ 
ical issue is then a reliable estimate of the uncertainty 
from A/;;;/, corrections, which remains a difficult and 
partly speculative problem. I refer to ll 103 . 1041 on this 
point and for many references on this very active field 
of LHCb physics. 

5. Summary and Outlook 

Factorization of exclusive and semi-inclusive B de¬ 
cays has developed into a mature theory, making a large 
amount of previously intractable final states accessi¬ 
ble to theoretically solid interpretations. The presence 
of energetic particles or jets implies many similarities 
with factorization methods applied in high-energy col¬ 
lider physics. A conceptually interesting point is that 


for exclusive decays the power-suppressed soft interac¬ 
tion = q s Wcil/)_ Lc g + h.c. related to the emission (or 

rather absorption) of soft quarks q c —> gi, c +q s appears at 
leading order in the expansion of the hard scale, whereas 
in collider physics only the familiar eikonal interaction 
gin- -Aj^-g describing soft gluon emission q c —> q c +g s 
is ever required at leading power. This is because the 
decaying B meson provides a “source” of soft partons. 

Beyond the leading power in the A /nib, the available 
theoretical tools often cease to be effective. For jet-like, 
semi-inclusive final states, it appears that soft-collinear 
factorization can in principle be extended to any order 
as long as the invariant mass of the inclusive hadronic fi¬ 
nal state is 0{mt,K). However, in practice already at the 
first sub-leading order in A /nib a large number of non- 
perturbative soft functions appears, as was discussed for 
the case of B —» X u tv. There is presently no method that 
could compute these functions, since lattice QCD can¬ 
not be used for the matrix elements of non-local opera¬ 
tors with exactly light-like field separations. For exclu¬ 
sive decays, a theory of power-suppressed effects sim¬ 
ply does not exist. 

Over the past ten years almost all calculations nec¬ 
essary to push the accuracy of charmless hadronic two- 
body decays to 0(a 2 ) have been performed, except for 
the two-loop QCD penguin amplitude. Once this is 
completed, direct CP asymmetries can for the first time 
be predicted with reliable perturbative uncertainties. In 
view of the upcoming BELLE II experiment an update 
of theoretical predictions with NNLO uncertainty and 
improved hadronic input parameters is timely. There 
are still many unmeasured, but theoretically predicted 
final states. Meanwhile, the LHCb experiment yields 
data on the electroweak penguin decays B —» 
the precision of which provides or will soon provide a 
challenge to theory. Extending the theoretical calcula¬ 
tions to 0(a 2 ) is difficult here, but since the spectator¬ 
scattering contribution is presently known only at its 
first order, the two-loop calculation should be done. 
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